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Motivation: Sign problem, Silver Blaze problem 
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Lefschetz-thimble path integral 


Finite-density QCD? 


Neutron star 
e Cold and dense nuclear matters 
© 2m, neutron star (2010) 
ə Gravitational-wave observations (2016~) 


Reliable theoretical approach to equation of state must be developed! 


Z(T,u) = f DA Det(P(A, u) +m) exp —Sym(A). 


quark gluon 
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Sign problem in finite-density QCD 
QCD & |QCD| 


ZQCD = [ra (det DP) e SYM | Ziacp| = [ea [det P| eM, 


If these two were sufficiently similar, we have no practical problems. 
However, it was observed in lattice QCD simulation that 
(e.g., Barbour et. al. (PRD 56 (1998) 7063)) 


ng 


" 
/ 
i 
ri 


-m4/2 -mn/3 


> Lg 
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Baryon Silver Blaze problem 


The curious incident of the dog in the night-time (Holmes, Silver Blaze). 
Problem: It is of great importance for finite-density lattice QCD 
to understand ng = 0 for fg < my/3. (Cohen, PRL 91 (2003) 222001) 


Current situation: For jig < m„/2, the problem is almost solved: 
Quark det. is (Adams, PRD 70 (2004) 045002) 


Det(P(A, ug) +m) | 
Det(D(A,0) +m) — 


I o60,-24. 


Re(AA) «uq 


and ess-min4(Re(A4)) = Mx/2 (Nagata et. al., PTEP 2012 014103). 


For fg > m, /2, no one knows how to understand this. 
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Method: Path integral on Lefschetz thimbles 
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Lefschetz-thimble path integral 


Complexification of fields 


There are two “new” approaches to the sign problem: 
e Complex Langevin method: Solve the Langevin eq. 


1] is a real stochastic noise, (n(0)n(0")) = 26(0 — 6’). 
e Path integral on Lefschetz thimbles: 


/ d'ea 99 — V (K,,R) / dag Pm 


oO 


Path integral is performed on steepest descent paths Js. 
This is mathematically rigorous, and always true for finite 
dimensional integrals. 
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Lefschetz thimble for Airy integral 
Airy integral is given as 


3 
Ai(a) = [ri (= + 2 


Complexify the integration variable: z = x + iy. 


Im(6) 


Integrand on R, and on A 
(a — 1) 
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Rewrite the Airy integral 
There exists two Lefschetz thimbles J, (c = 1,2) for the Airy 


integral: 
Ai(a) Son / “an am 
= 5 — — az]. 
= J, 2T p 3 


ng: intersection number of the steepest ascent contour Ke and R. 
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Lefschetz decomposition formula 


Oscillatory integrals with many variables can be evaluated using the 
“steepest descent” cycles 7,: (classical eom S’(z,) = 0) 


/ d"z e 3) = > (Ke; 3 d'zg 9i 


o 


Js are called Lefschetz thimbles, and Im[S] is constant on it: 


M =}, dzi(t) _ eo 


1—-—oo dt pe 


= [40 


(Ko, R}: intersection numbers of duals K, and R” 


(Ko = {2(0)|2(00) = zo }). 


[Witten, arXiv:1001.2933, 1009.6032] 
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Analysis: Semi-classical analysis of the one-site Hubbard model 
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Lefschetz-thimble path integral 


One-site Fermi Hubbard model 
One-site Hubbard model: 


A = Uff, — u(fu + f). 


Fock state gives the number density immediately: 


"d 


(YT, Hidaka, Hayata, 1509.07146)(cf. Monte Carlo with 1-thimble approx. gives a wrong result: 


Fujii, Kamata, Kikukawa,1509.08176, 1509.09141; Alexandru, Basar, Bedaque,1510.03258.) 
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Analysis 


Path integral for one-site model 


The path-integral expression for the one-site Hubbard model: : 


B n A 2. UD 
Z24 — | ay (14 efe)? 78/20 
2nU R mar). 


Fermion Det 


Integrand has complex phases causing the sign problem. 


ọ is an auxiliary field for the number density: 
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Analysis 


Flows at u/U < —0.5 (and u/U > 1/5) 


0.75 


0.25 


0 
Re(z) 


Figure: Flow at u/U=-1. Ix 


-0.25 


-0.5 


cR. 


N dz e- 909, 


Number density: n, = 0 for u/U < —0.5, n, = 2 for u/U > 1.5. 


(YT, Hidaka, Hayata, 1509.07146) 


Z- 
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Analysis 


Flows at —0.5 < u/U < 1.5 
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Figure: Flow at u/U = 0 


-0.75 


py p 1. 


This value is far away from n = Im (z)/U=0, 1, or 2. 


m 


Complex saddle points lie on Im(z,,)/U 
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Complex classical solutions 


If GU > 1, the classical sol. 
for —0.5 < u/U < 1.5 
are labeled by m € Z: 


U 
Ea cR By +2rmT. 


At these solutions, the classical actions become 


Sox -£ (+, 
Re (Sin — So) ~ 2m m?, 
Im S. c 2mm (5 4- 1). 


Yuya Tanizaki (RIKEN, RBRC) Lefschetz-thimble path integral May 26, 2016 @ RBRC 


u[U x 1/2 


17 / 36 


Semiclassical partition function 


Using complex classical solutions Zm, let us calculate 


oo 
Za t= > p nm. 


M=—-09 


This expression is valid for —1/2 S u/U S 3/2. 
This is calculable using the elliptic theta function: 


= 1 
-S [149 9 Heo + \a—27?m?/BU 
e ( + ) cos 27m (5 + 2 e 


Zel = 
m=1 
a H 1 2 2 
= gr |r| 4 |,6 277" }. 
U 2 
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Number density & Lee-Yang zeros 
Lee-Yang zeros of Za: 


e o% 
. è 
r . Re 
° e 


Semiclassical study gives the correct transition! 


2 (1< u/U «2. 
InZa4—4 1 (0<p/U « 1), 
0 (—1/2 « u/U « 0). 


(YT, Hidaka, Hayata, arXiv:1509.07146[hep-th]) 
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Important interference among multiple thimbles 


Let us consider a “phase-quenched” multi-thimble approximation: 


Za) = Do le | = e 9 45(0, e 0727). 


e Lee-Yang zeros cannot appear at u/U = 0, 1. 
e One-thimble, or “phase-quenched”, result: n ~ u/U + 1/2. 


Consequence 


In order to describe the step functions, we need interference of 
complex phases among different Lefschetz thimbles. 


(cf. Particle Productions: Dumulu, Dunne, PRL 104 250402) 
(cf. Hidden Topological Angles: Behtash, Sulejmanpasic, Schafer, Unsal, PRL 115 041601) 
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Numerical results 


Results for GU = 30: (1, 3, 5-thimble approx.: Jo, Jo U Ji, and Jo U J41 U Ja ) 


2.0 
— exact 
— 1-thimble 


15 --- 3-thimble 
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Necessary number of Lefschetz thimbles ~ 8U/ (27). 


(YT, Hidaka, Hayata, NJP 18 (2016) 033002, arXiv:1509.07146[hep-th]) 
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Discussion 


Discussion: Speculation on Silver Blaze problem 
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Discussion 


Curious incident of n in the one-site model 


We have a big difference bet. the exact result and naive expectation: 


05 O 0.5 1 15 ajo 


Naive guess: This is similar to what happens for QCD and |QCD|. 
p/U = -0.5 & u =m, /2. 
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Dimensionally reduced expression of fermion det. 


One-site Hubbard model: 


U 
Det lo — (» + 5 + 3] = (1 a e &CUA- nolle? l 


Quark determinant in QCD: 


Det aD, +m) — p] = N (A) T] terem) Lee Perret), 


£j20 
where the spectrum of yu(D, + m) is 
Aga) = &j( A) — i$; (A) + (2n + 1)izT. 


Formal correspondence: —U/2 + ess-min(e;) = m4/2, o € $;. 
(YT, Hidaka, Hayata, NJP 18 (2016) 033002, arXiv:1509.07146[hep-th]) 
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Discussion 


Silver Blaze problem for u < —U/2, u < m4/2 


One-site Hubbard model: As BU >> 1 and —U/2 — u > 0, 


U 
Det lo = (» + 7 + e) = (1 + EURER)" ~]. 


The sign problem almost disappears, so that J, ~ R. 
Finite-density QCD: As 6 — oo and u < m,/2, 


Det ud, +m) — ul 
Det l4 D 4 F m)| 


The sign problem disappears by the reweighting method. 
=> Lefschetz thimbles ~ Original integration regions 


> 1 
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Silver Blaze problem for u > —U/2, u > m,/2 


One-site Hubbard model: At each real config., the magnitude is 
exponentially large: 


Det lo — (» + 5 + e) =.) (eP U +e/2)) 


This large contributions must be canceled exactly in order for n = 0. 


In the Lefschetz-thimble method, this is realized by interference 
among multiple thimbles. 


Question 
Does the same thing happens for QCD? | 


We guess that it would, but we do not yet know how severe it is. 
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Bonus 


Bonus: Complex Langevin method 
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Complex Langevin method 


Complex Langevin has been regarded as a sign-problem solver via 
stochastic quantization (Klauder, PRA 29, 2036 (1984), Parisi, PLB 131, 393 (1983)): 


dz(0) _ 8S 
0: Stochastic time, 7: Random force (n(0)n(0')),, = 20(0 — 0"). 
Ito calculus shows that 


x (Olen(0)) = A(O" (z,(8))), — (0 (24(8)) S (27(0)))n: 


If the I.h.s becomes zero as 0 — oo, this is the Dyson-Schwinger eq. 
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Bonus 


Complex Langevin and Lefschetz thimbles 


For any solutions of the DS eq, 


(Olh = 2 374. | ase oq), 


c o 


ind, € C. To reproduce physics, d, = (Ko, R) € Z. 


So far, we ONLY assume the convergence of the complex Langevin 
method. 
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Bonus 


Semiclassical limit 
Let us take h < 1 for computing 


(Oln) = ds [ dzeS/*O(z). 


| have NO idea how to compute the LHS. However, positivity of the 
probability density and its localization around z's imply that 


"ce. 3 (Ole) S 2.508). 


RHS is 


[ dze SOM"O(z) > S Ole). 
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Bonus 


Semiclassical inconsistency 
In the semiclassical analysis, one obtains (for dominant saddle points) 


m NNNM 
Z \ Be.) l 


Co 


Co > 0. And, d, = (Ks, R) to get physics. = Inconsistent! 
(Hayata, Hidaka, YT, arXiv:1511.02437[hep-lat]) 
We show that the complex Langevin is wrong if 
e There exist several dominantly contributing saddle points. 
e Those saddle points have different complex phases. 


Open question: Recently, our knowledge on the CL method is 
developing (Aarts, Seiler, and Stamatesc, 2009; Nishimura, Shimasaki, 2015, etc.). 
Can one understand our result from this viewpoint? 
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Bonus 


Proposal for modification 
Assume that 


_ Ko, R) 2nh e- S62. 


S" (x. 


Because of the localization of probability distribution P, it would be 
given as 


P= GPs supp(P,)Nsupp(P;) = 0. 


E 


Assumption means “CL = phase quenched multi-thimble approx.” : 


(Ko, R) 


2h 
T e SGo)/h 


Sz) De): 


(0(2)) 2», 7 


o 
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Bonus 


Proposal for modification (conti.) 


If so, defining the phase function 


|S" (ze)| —iIm S(z;,)/h = 
S" (Ze) E id Xsupp(Ps) (2, Z), 


we can compute 


(©(2,, Zn)O(Zn))n 
((25, 24))q l 


This new one is now consistent within the semiclassical analysis. 
(Hayata, Hidaka, YT, arXiv:1511.02437[hep-lat]) 

Caution: Our proposal evades inconsistency, but is not necessarily 
correct. Can we improve the proposal? 


(0()y = 
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Bonus 


Complex Langevin study of one-site Hubbard model 


BU =30 u/U=0 


o 


o 
bo 


8[Im(z) — (u + U/2)] 
S 


= 
io 


=S  -10  -5 0 5 10 15 
BRe(z) 
Figure: CL distribution at u/U = 0 


It looks quite similar to Lefschetz thimbles. 
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Bonus 


Complex Langevin simulation 
One-site Fermi Hubbard model: 


=== exact 
2.0} HEH w/o phase 


‚9: w/ phase 

= 1.5 
o 
Ei — — — | |$*X95Emut5. .. 

0.5 

0.0 

-0.5 0.0 0.5 10 Ls 
p/U 


(Hayata, Hidaka, YT, arXiv:1511.02437[hep-lat]) 


Consequence 


Modified complex Langevin is not perfect yet, but it seems to point a 
correct way. 
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Summary & Conclusion 


Summary and Conclusion 


e Picard-Lefschetz theory gives a suitable framework for 
saddle-point analysis even if S(@) takes complex values. 


e One-site Hubbard model is a nice toy model to play with the 
sign problem. 


e Destructive and constructive interference of complex phases 
among Lefschetz thimbles play a pivotal role for the baryon 
Silver Blaze. 
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